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Abstract. In this talk I address three topics related to the shape of hadrons: 

1. The Lorentz contraction of bound states. Few dedicated studies of this exist - 1 describe a recent calculation for ordinary 
atoms (positronium). 

2. Does the A-dependence of nuclear structure functions indicate a change of proton shape in the nuclear environment? 
(My short answer is no.) 

3. The size of Fock states contributing to processes involving large momentum transfers. End-point configurations can be 
transversally extended and yet sufficiently short-lived to contribute coherently to hard scattering. 

LORENTZ CONTRACTION OF BOUND STATES 

We are all familiar with the classical prediction of special relativity that rigid rods appear contracted to a moving 
observer. Lorentz contraction is also often appealed to in descriptions of the scattering of hadrons and nuclei, which 
are depicted as thin disks when moving with high momenta. It is thus surprising that few studies seem to have been 
devoted to the contraction of bound state wave functions in relativistic quantum physics. I am aware of only a few 
examples [1], most of which are in the context of 1+1 dimensional models. Why do we not, in courses on relativistic 
quantum mechanics, derive the contracted wave function of the ordinary hydrogen atom in relativistic CM motion? 
What do we in fact know about the Fock states of the atomic wave function (even at lowest order in a) when the 
(kinetic) energy of the atom is much larger than the masses of its constituents? Here I shall briefly describe some 
results of a recent study by my student Matti Jarvinen [2] of these issues, which are quite fascinating and deserve more 
attention. 



Wave functions and Fock states 

The wave function of a bound state can be defined in several ways, depending on how the time slice (snap-shot) 
is taken of the propagating state. In non-relativistic quantum mechanics we are used to equal time wave functions 
1Po(*; x i) . . . ,x^), which are probability amplitudes for finding the N constituents at positions Xi, . . . ,X# at a common 
time t (see Fig. la). The concept of equal time is frame dependent, which makes boosting ^Po wave functions non- 
trivial. 



1 Invited talk at the Workshop on the Shape of Hadrons, 27-29 April 2006 in Athens, Greece. 
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FIGURE 1. (a) The equal time wave function is the amplitude for finding the constituents of a bound state at relative distance 
x when measured at a common time t. (b) The light-front wave function x ¥ + is the amplitude for finding the constituents of a bound 
state at relative distance (x~ , ) when measured at a common LF time x + = t + z- 



In relativistic descriptions wave functions v F+(jt + ;jt]~,xix, . . . ,xZ,Xnx) defined at equal light front (LF) time 
x + = t + z are commonly used [3] (Fig. lb). Each constituent is then observed at a common x + = x\ = . . . = xt, 
and the x ¥ + wave function describes their distribution in x~ =t — z and x± = (x,y). This LF wave function is invariant 
under boosts along the z-axis. Such boosts leave the (relative) LF time invariant: x' + = e^x + , so that x + = transforms 
into x /+ = for any boost parameter £. 

When the constituents move non-relativistically it makes no difference whether they are observed at equal t or equal 
x + : 'Pq = x ¥ + since the constituents move only an infinitesimal distance in the time it takes a ray of light to connect 
them (Fig. lb). However, when an internally non-relativistic system like the Hydrogen atom is in relativistic motion 
also its constituents move relativistically. Then the two definitions of the wave function are no longer equivalent. 
Which one should be expected to exhibit Lorentz contraction? 
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FIGURE 2. Lorentz contraction for the equal time wave function Wo[t,x), as it is commonly depicted. 

In the classical derivation of Lorentz contraction each observer measures the length of the rod when its ends are at 
the same (ordinary) time t. Contraction is the result of the different perception of time for observers in relative motion. 
It is thus clear that we should expect contraction for the equal (ordinary) time wave function (Fig- 2). As we already 
noted, V P+ is in fact invariant under boosts. 

The pictorial representation of contraction in Fig. 2 leaves, however, many open questions. A field theory bound 
state is defined as the solution of the Schrodinger equation 



H\N(P,s)) =E(P)\N(P,s)) 
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where the hamiltonian H contains operators that create and annihilate particles 2 . The bound state is a superposition of 
Fock states, e.g., for the proton 
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Should all the various Fock components contract at a universal rate? This would imply a simple relation between equal- 
time wave functions in different frames - which is not expected since boosts are dynamical operators [4]. Boosting 
a rest frame wave function is in fact as complicated as solving the Schrodinger equation (1) directly in the boosted 
frame. 



- Since we are dealing with motion where the kinetic energies exceed particle masses we cannot restrict ourselves to Fock states with a fixed particle 
number. Even for the Hydrogen atom, Fock states beyond e + e~ can be neglected only in the rest frame (and at lowest order in a). 
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FIGURE 3. Momenta which are equal and opposite in the rest frame can be boosted into the forward hemisphere. 



A Lorentz contraction changes the magnitudes of longitudinal momenta, leaving the transverse momenta un- 
changed. In particular, contractions do not change the sign of longitudinal momenta. Boosts on the other hand can 
readily turn backward moving particles into forward moving ones, as indicated in Fig. 3. Thus it would be miraculous 
if bound state wave functions only would contract under boosts. 
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FIGURE 4. In old-fashioned perturbation theory there is a separate diagram for each time ordering of the interaction vertices. 
Only the sum of all diagrams is Lorentz invariant. 

It is also well-known that Fock states mix under boosts. A simple 3 example is shown in Fig. 4 using time-ordered 
(old-fashioned) perturbation theory. The two diagrams differ by the time ordering of their interaction vertices. The first 
diagram of Fig. 4 has a single particle in the intermediate state, while the second one has five. Their relative magnitude 
is frame dependent, but their sum is Lorentz invariant, 
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where Ej = y p? +m 2 and E — \J (pi +P2) 2 +M 2 . In the limit of very large boosts (i.e., in the infinite-momentum 
frame) only the first diagram survives - and the rules are then the same as those of LF time ordered perturbation theory. 



Positronium in relativistic motion 

A relativistic bound state has (infinitely) many Fock components, making it practically impossible to derive the 
equal time wave function in the rest system, and even less so in a general frame. However, the Hydrogen atom (or 
positronium, for simplicity) has a very simple wave function in the rest frame. At lowest order, only the \e + e~) Fock 
state of positronium contributes, bound by the instantaneous Coulomb potential. The internal motion is non-relativistic, 
suggesting essential simplifications also in a frame where the bound state moves with relativistic speed. Thus we may 
ask whether the wave function of the \e + e~) component is Lorentz contracted? Do higher Fock states contribute, and 
do their wave functions also contract? 

I shall here only summarize the results obtained in [2]. It is necessary to start from the fully relativistic framework 
provided by the Bethe-Salpeter (BS) equation (Fig. 5). By inspecting the equation one may surmise that many 
contributions, such as electron pair production and annihilation, are of higher order in a and thus can be neglected in all 
frames. The relevant terms remaining in the interaction kernel are instantaneous photon exchange (as in the rest frame) 
and transverse photon exchange. Thus at lowest order the BS equation simplifies to that shown in Fig. 6. As might have 
been expected, the boost has given the Coulomb exchange of the rest frame a transverse, propagating component. The 
propagation time is of the order of the Bohr radius, 1 /(am e ). Since the binding energy of positronium AE « a 2 m e , the 
lifetime of the \e + e~) Fock state is longer by a factor 1 /a compared to the lifetime of the |<? + e~y) state. The likelihood 



FIGURE 5. The Bethe-Salpeter equation. The blobs represent the positronium wave function 'Fp, K is the interaction kernel and 
S is the e + e~ propagator. 

of "seeing" a photon in flight is thus of G (a), even though transverse photon exchange contributes to positronium 
binding at leading order. The probability of catching two photons in flight is of G (oc 2 ) and can be neglected. 
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FIGURE 6. The time-ordered bound state equation for positronium in the ladder approximation. The blobs denote the equal time 
wave function <pp of the e + e~ Fock component carrying total 3-momentum P. 

The leading order BS equation shown in Fig. 6 can be readily solved. Due to the time ordering explicit Lorentz 
symmetry is lost. Hence there is no guarantee that the positronium energy is correctly given as E = VP 2 + M 2 , where 
M is the rest mass - except that we know that QED is Poincare invariant. Indeed, this dependence of the energy on the 
CM momentum P emerges almost miraculously from the equations. The \e + e~) Fock amplitude <jjp satisfies 
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where Ej, = VE 2 — P 2 — 2m e , y = y^P 2 + (2m e ) 2 /2m e and p is the relative momentum between the electron and the 
positron. The scaling pu — > ypn and h\ — > yk\\ eliminates all explicit P-dependence from the bound state equation. This 
means that Ej, is independent of P as it should be, and that <pp Lorentz contracts. 

The wave function of the \e + e~ y) Fock component can be similarly found. The probability distribution of the photon 
(integrated over the electron momenta) turns out to factorize for the positronium ground state, 
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where k t = (kj_,fci/y) is the photon momentum corrected for contraction and 9 C is the angle between P and k c . The 
distributions are 
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FIGURE 7. The angular dependence of the contracted and integrated photon distribution (5) in the positronium ground state. The 
lines show the angular distribution /(cos 0)/(yj3 2 ) [defined in (6)] for j8 = 0.001, 0.5, 0.9 and 0.999. For j3 = 0.001 (solid line) 
the distribution is close to the symmetric limit of the rest frame. For /3 = 0.999 (dotted line) the distribution approaches the infinite 
momentum frame (alias LF wave function) where the momenta of all particles are in the forward hemisphere. 



where <f>o is the \e + e~) Fock amplitude in the rest frame (P = 0), and j3 = — 1 /y 2 . The angular distribution (6) of 
the photon depends on the CM motion (j3 , y) even after rescaling. Thus this Fock component does not simply Lorentz 
contract, as we already anticipated (cf. Fig. 3). As shown in Fig. 7, the photon distribution changes with increasing /3 
from being symmetric (for j5 ~ 0) into one that for j8 — s- 1 has support only in the forward hemisphere (cos 9 C > 0). The 
latter (infinite momentum frame) distribution agrees with that of the wave function defined at equal LF time x + =t + z 
(which is invariant under boosts). 

The study [2] of the frame dependence of the Hydrogen (or positronium) equal time wave function has thus revealed 
a richness beyond the standard picture of Fig. 2. The wave function of the lowest \e + e~) Fock state (which is non- 
relativistic in the rest frame) indeed Lorentz contracts. Note, however, that this contraction is in the relative longitudinal 
momentum of the constituents. In a positronium state of momentum P both the electron and the positron momenta 
equal P/2, up to the small relative momentum of & (a). For P = (the left hand side of Fig. 2) the momenta of the 
constituents are in opposite hemispheres. For P larger than the Bohr momentum both the electron and the positron 
momenta are in the forward hemisphere, the flattened disk on the right hand side of Fig. 2 describing only their 
deviations from P/2. 

For the photon in the \e + e~y) Fock state the situation is different since it carries only a fraction (a) of the bound 
state momentum P. As seen in Fig. 7, the photon momentum can be in the backward hemisphere even for large P, 
being restricted to the forward hemisphere only in the infinite momentum frame (P — > °°). 

Further information on frame dependence could be obtained from studies of the higher order corrections to the 
positronium wave function including, e.g., Fock states with an additional fermion pair. It should also be possible to 
study the CM momentum dependence of non-relativistic bound states in theories with other (e.g., 3 ) interactions. The 
frame dependence of a state such as the proton, which is relativistic also in its rest frame, is likely to be more involved 
than that of positronium. It is not obvious that any Fock component of the proton (including the lowest one, \uud)) 
obeys the classical longitudinal contraction sketched in Fig. 2. 



SHAPE OF NUCLEONS IN NUCLEI 

The discovery of a nuclear (A-)dependence of parton distributions raised the question whether the cause is to be found 
in an actual modification of the nucleon wave function in the nuclear environment, or reflects the dynamics of the 
Deep Inelastic Scattering (DIS) measurement [5]. To some extent this is a matter of definition. The parton distribution 




FIGURE 8. The partem distribution (8). Rescattering arising from the Wilson line (9) is shown as dashed lines. 



shown in Fig. 8, 
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(where x =t±z) involves the Wilson line 



W[x~,Q] = Pexp 
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which may be included in the definition of the LF wave function (which then becomes specific to DIS and similar hard 
processes). On the other hand, the Wilson line may also be regarded as reflecting physical rescattering of the struck 
quark on target spectators [6] (dashed lines in Fig. 8). Since there are more spectators in a nuclear environment the 
rescattering can be expected to be A-dependent. This view of rescattering is the same as the one usually formulated in 
the "target rest frame" in models of shadowing [7, 8]. In either case, the rescattering must be coherent with the hard 
virtual photon interaction in order to affect the total DIS cross section. This coherence allows the different physical 
interpretations to be equivalent. 
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FIGURE 9. Schematic behaviour of the ratio of the nuclear to free nucleon Fi structure functions R^(x B ,Q 2 ) as function of 
Bjorken x B for a given fixed Q 2 . From Ref. [8]. 

As sketched in Fig. 9, the A-dependence extends beyond the shadowing region of small x B . One may ask whether the 
"antishadowing" and "EMC" effects reflect an actual change of the nucleon shape, beyond effects due to rescattering 
of the DIS probe. One answer is provided by the observation made several years ago [9] that these structures essentially 
vanish if one views the parton distribution in coordinate (xT) space, rather than in momentum (x B ) space. The parton 
distribution (8) is a Fourier transform of an overlap of nucleon wave functions, with a shift x~ in the position of the 
struck quark (Fig. 10). We may express the the matrix element in terms of the parton distributions through the inverse 
Fourier transform, 

(N\q(x')Y + W[x^ ,0}q(0)\N) - (x~ -> -x~) = Aim J dx B [f g / N {x B ) + fg/ N (x B )] sin (\mx B x^) (10) 




FIGURE 10. The matrix element in the integrand of the partem distribution (8), i.e., the handbag diagram of Fig. 8 viewed in 
coordinate space (rescattering is not shown). The position of the struck quark differs by x~ in the two wave functions (whereas 
x+ = %j_ = 0). 



The rhs. of this equation is essentially given by the Fi structure function. Thus we can study the A -dependence of the 
parton distribution in coordinate space, defined as 

cf{x-,Q 2 )^ [ l ^F°(x Bl Q 2 )R A F2 (x Bl Q 2 ) S m(±mx B x-) (11) 
Jo x B 

where R F (xb, Q 2 ) is the experimentally measured ratio of nuclear to deuterium structure functions sketched in Fig. 9. 
The corresponding ratio in coordinate space, defined as 

r ^^9(^¥) (12) 

can then be formed using data on structure functions and is shown in Fig. 1 la. 
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FIGURE 11. (a) The coordinate space ratio R A (x~ ,Q 2 ) (12) obtained by Fourier transforming data on F^ixBiQ 2 ) structure 
functions for A = He, C and Ca. (b) The momentum space ratio R c (xb,w, Q 2 = 5 GeV 2 ) for Carbon, obtained by Fourier 
transforming a modified coordinate space distribution in which all nuclear effects are eliminated for x~ < w. 

Within the ca. 1% error bars [9] the ratio R A (x~ ,Q 2 ) is consistent with having no A-dependence for x~ < 5 fm. At 
longer distances x~ > 5 (i.e., t = — z > 2.5 fm since x + — 0) shadowing sets in. Thus viewed from coordinate space 



(which is not unnatural for discussing effects of nuclear size) we may regard 3 antishadowing and the EMC effect as 
merely resulting from Fourier transforming a fiat distribution (of finite length) in x~~ ! This is corroborated in Fig. 1 lb, 
where the reverse transform back to momentum {xb-) space is made, under the assumption that R A (x~ ,Q 2 ) is unity 
for x~ < w (and takes the values of Fig. 11a for x~ > w). It is seen that the antishadowing and (most of) the EMC 
effect is reproduced assuming no nuclear dependence in coordinate space for x~ < 5 fm. The nuclear effects can thus 
be ascribed solely to shadowing. 

The parton distribution ^(x - ,Q 2 ) in coordinate space is insensitive to the region of Fermi motion at large Xb in 
Fig. 9, where the structure function F2{xb 1 Q 2 ) is small. The sizeable nuclear dependence of R^- 2 (xb,Q 2 ) at large xb 
reflects the ratio of very small which do not appreciably affect the inverse Fourier transform (11). 



SIZE OF HARD SUBPROCESSES 

The third aspect of shape that I would like to discuss concerns the size of coherent hard subprocesses in scattering 
involving large momentum transfers. As sketched in Fig. 12, in inclusive DIS {ep — > eX) we expect that the virtual 
photon (whose transverse coherence length is ~ 1/2) scatters off a single quark. The quark is typically part of a Fock 
state with a hadronic, ~ 1 fm size. In elastic scattering {ep — > ep), where the entire Fock state must coherently absorb 
the momentum, one might on the other hand expect [11] that only compact Fock states of the photon, with transverse 
sizes r± ~ l/Q will contribute. Thus the dynamics of inclusive and exclusive processes appears to be quite different. 
In particular, the dependence on the electric charges of the quarks is expected to be, qualitatively, 
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FIGURE 12. The virtual photon scatters from single quarks in inclusive deep inelastic scattering {left). If the valence quarks 
absorb equal shares of the momentum transfer in the exclusive ep — > ep process (right) only compact Fock states can contribute. 

In contrast to these expectations the data suggests a close connection between inclusive and exclusive scattering. 
The resonance production ep — * eN* cross sections (including N* = p) average the DIS scaling curve when plotted at 
the same value of xb (or of the related Nachtmann variable |) [12]. Examples of this Bloom-Gilman duality are shown 
in Fig. 13. A natural explanation of duality is that the same Fock states of the proton contribute in both cases [13]. 
Resonance formation occurs on a longer time scale than the hard subprocess, hence is incoherent with it and cannot 
change the total cross section. Only the local mass distribution (resonance bumps) is sensitive to the hadronization 
time scale. 



3 Understanding the dynamics of nuclear dependence in momentum space is nevertheless interesting in its own right. See [10] for recent ideas about 
the origin of the antishadowing enhancement. 




FIGURE 13. Bloom-Gilman duality, (a) Resonance production in exclusive ep — > eN* average the scaling curve measured at 
high Q 2 as a function of the Nachtmann variable % (i.e.,xg with a target mass correction), (b) Proton Fi structure function in the A 
(top) and Si i (bottom) resonance regions from Jefferson Lab Hall C, compared to the scaling curve from Ref. [16]. The resonances 
move to higher § with increasing Q 2 , which ranges from ~ 0.5 GeV 2 (smallest £ values) to ~ 4.5 GeV 2 (largest £ values). From 
Ref. [12]. 



As may be seen from the rhs. of Fig. 13 the resonance bumps migrate with increasing Q 2 to large xg, where the 
struck quark carries most of the proton momentum. This is the "end-point" region where the Lepage-Brodsky picture 
[11] of exclusive form factors breaks down [14]. When the struck quark carries nearly all the proton momentum the 
remaining (soft) partons can join the wave function of the final proton without the hard gluon exchanges indicated on 
the rhs. of Fig. 12 (this is the 'Feynman mechanism'). The soft partons have short formation times <*; 1 — xg and can 
thus stay coherent with the hard process, even though they are broadly distributed (r± ~ 1 fm) [15]. 

The correct dynamics of exclusive form factors is still an open and important issue. Further hints may be obtained 
from spin dependence. Here I would like to mention the photoproduction process yp — ► pp, which has been measured 4 
at large momentum transfer \t\ < 10 GeV 2 at HERA [17, 18]. Due to the high energy one expects two-gluon exchange 
to dominate. The upper vertex (Fig. 14) should be hard due to the large momentum transfer \t \ carried by the gluons. 
Consequently the quark helicity is conserved at the gluon and photon vertices and the sum of the quark and antiquark 
helicities vanish, as indicated in Fig. 14. In the absence of orbital angular momentum (which should be suppressed 
due to the small transverse size of the quark pair) the p meson will then be longitudinal (X p = 0). 
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FIGURE 14. Upper vertex of the photoproduction process ep —* pp. The photon creates a quark and an antiquark which carry 
fractions z and 1 — z of its momentum. The pair scatters with large momentum transfer via two-gluon exchange and then forms the 
p meson. The transverse size r ± of the quark pair is of ^(1/ yjfj) unless z —* or z — > 1. 




4 The final proton was not observed and may dissociate into a low mass system. 



The data shows that the p is dominantly produced with helicity X p = ±1. The density matrix element 
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where the Tx p x y are helicity amplitudes, is small in the whole range of \t\ ~ 6 GeV (Fig. 15). Thus quark helicity is 
maximally violated, while the y and p helicities are approximately equal. 
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FIGURE 15. The spin density matrix element rSn of (14) for p meson photoproduction as a function of Itl from HI [18] (full 
points) and ZEUS [17] (open points). The dotted line shows the prediction of a two-gluon exchange model. The dashed line shows 
the expectation from s-channel helicity conservation (SCHC). From Ref. [18]. 

The perturbative y+ (gg) — > qq with a quark helicity flip such that the helicity of the quark pair is X qq = 1 is [19] 
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This amplitude should be convoluted with the distribution amplitude of the p meson, which typically has the form 

4>p(z)c<z(l- z ) (16) 

The integral over z is then logarithmically singular at z = 0,1, indicating the potential importance of endpoint 
contributions [19]. 

When a quark carries low longitudinal momentum its transverse size distribution grows rapidly (due to the absence 
of Lorentz time dilation). It may then scatter softly and flip its helicity without the m q /\f^-i penalty of (15). This 
dynamics is actually quite similar to that of ordinary DIS viewed in the target rest frame. The virtual photon splits into 
a quark and an antiquark, with the antiquark carrying finite momentum in the target rest frame (even though the photon 
energy v — > °°). The soft scattering of the antiquark in the target "releases" the quark, which forms the current jet. In 
the photoproduction process yp — > pp the distribution amplitude (16) of the p meson suppresses the endpoint region, 
but apparently not sufficiently. Thus the dominant subprocess may actually have transverse size r± ~ 1 fm rather than 
r ± -1/7^7. 
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